expre<-read.table("out_sig_eqtl_gene_expr",sep="\t")
geno<-read.table("out_sig_eqtl_genotyping",sep="\t")
eqtl<-read.table("out_only_sig_eqtl",sep="\t")
pdf("signi_eqtl_gene_expr.pdf")
hist(expre[,4],breaks=seq(0,0.05,0.0001),xlim=c(0,0.05),xlab="Pvalue",ylab ="counts",main="signi_eqtl_gene_expr_survival_pvalues")
dev.off()
pdf("sig_eqtl_genotyping.pdf")
hist(geno[,4],breaks=seq(0,0.05,0.0001),xlim=c(0,0.05),xlab="Pvalue",ylab ="counts",main="sig_eqtl_genotyping_survival_pvalues")
dev.off()
library(gplots)
pdf("venn1_snps.pdf")
list1<-list(unique(geno[,2]),unique(expre[,2]))
names(list1)<-c("sign_gene_surv_eqtl_snps","sign_snps_surv_eqtl_snps")
venn(list1)
dev.off()
pdf("venn1_gene.pdf")
list2<-list(unique(geno[,1]),unique(expre[,1]))
names(list2)<-c("sign_gene_surv_eqtl_genes","sign_snps_surv_eqtl_genes")
venn(list2)
dev.off()
pdf("venn2_snp.pdf")
list3<-list(unique(eqtl[,2]),unique(geno[,2]))
names(list3)<-c("sign_eqtl_snps","sign_snps_surv_eqtl_snps")
venn(list3)
dev.off()
pdf("venn2_gene.pdf")
list4<-list(unique(eqtl[,1]),unique(expre[,1]))
names(list4) <- c("sign_eqtl_genes","sign_gene_surv_eqtl_genes")
venn(list4)
dev.off()
